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Abstract 



In this paper, detection of the primary user (PU) signal in an orthogonal frequency division multi- 
plexing (OFDM) based cognitive radio (CR) system is addressed. According to the prior knowledge of 
the PU signal known to the detector, three detection algorithms based on the Neyman-Pearson philosophy 
are proposed. In the first case, a Gaussian PU signal with completely known probability density function 
(PDF) except for its received power is considered. The frequency band that the PU signal resides is also 
assumed known. Detection is performed individually at each OFDM sub-carrier possibly interfered by the 
PU signal, and the results are then combined to form a final decision. In the second case, the sub-carriers 
that the PU signal resides are known. Observations from all possibly interfered sub-carriers are considered 
jointly to exploit the fact that the presence of a PU signal interferers all of them simultaneously. In the 
last case, it is assumed no PU signal prior knowledge is available. The detection is involved with a search 
of the interfered band. The proposed detector is able to detect an abrupt power change when tracing 
along the frequency axis. 
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I. Introduction 

Radio spectram is the medium for all types of wireless communications, such as cellular phones, 
satelHte-based services, wireless low-powered consumer devices, and so on. Since most of the usable 
spectrum has been allocated to existing services, the radio spectrum has become a precious and scarce 
resource, and there is an urgent concern about the availabiUty of spectrum for future needs. Nonetheless, 
the allocated radio spectrum today is not efficiently utilized. According to a report of the United States 
Federal Communications Commission (FCC) [1], there are large temporal and geographic variations in the 
utilization of allocated spectrum ranging from 15% to 85%. Moreover, according to Defense Advanced 
Research Projects Agency (DARPA), in the United States, only 2% of the spectrum is in use at any 
moment. It is then clear that the solution to the spectrum scarcity problem is dynamically looking for the 
spectrum "white spaces" and using them opportunistically. Cognitive radio (CR) technology, defined first 
by J. Mitola [2], [3], is thus advocated by FCC as a candidate for implementing opportunistic spectrum 
sharing. The spectrum management rule of CR is that all new users for the spectrum are secondary 
(cognitive) users and requiring that they must detect and avoid the primary user. 

To achieve the goal of CR, it is a fundamental requirement that the cognitive user (CU) performs 
spectrum sensing to detect the presence of the primary user (PU) signal. Digital signal processing 
techniques can be employed to promote the sensitivity of the PU signal sensing. Three commonly adopted 
methods are matched filtering, energy detection [4]-[10], and PU signal feature detection with the cyclo- 
stationary feature most widely adopted [11]-[14]. Moreover, cooperation among CUs in spectrum sensing 
can not only reduce the detection time and thus increase the agility, but also alleviate the problem that 
a CU fails to detect the PU signal because it is located at a weak-signal region [8]-[10], [15]-[20]. For 
overview of these approaches and their properties, see [21]-[23]. 

It is concluded in [24] that orthogonal frequency division multiplexing (OFDM) is the best physical 
layer candidate for a CR system since it allows easy generation of spectral signal waveforms that can fit 
into discontinuous and arbitrary-sized spectrum segments. Besides, OFDM is optimal from the viewpoint 
of capacity as it allows achieving the Shannon channel capacity in a fragmented spectrum. Owing to 
these reasons, in this paper, we conduct spectrum sensing in an OFDM based CR system. 

In detection theory, the Neyman-Pearson (NP) criterion is used when there is difficulty in determining 
the prior probabiUties and assigning costs [25] for hypotheses, which is the case in our PU signal detection. 
The NP detector compares the UkeUhood ratio (LR) with a threshold determined by the constraint of 
false alarm probability to decide which hypothesis is true. However, in many cases, some PU signal 
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parameters, such as power level, correlation properties, frequency band, and so on, may not be known. 
At this moment, the PU signal detection problem becomes a composite hypothesis testing, which requires 
performing estimation for those unknown parameters in the probability density function (PDF) of the 
observation for either hypothesis. Thus, the degree of detector complexity is directly related to the 
knowledge of the signal and noise characteristics in terms of their PDFs. Moreover, since the estimation 
error is not negligible, the detection performance decreases as we have less specific knowledge of the 
signal and noise characteristics. According to the prior knowledge about the PU signal, three cases of 
PU signal detection in a cognitive OFDM system are considered in this paper. 

In Case A, we assume the PU signal model is known and consider a Gaussian PU signal with completely 
known PDF except for its received power. The normalized, i.e. unity diagonal elements, covariance matrix 
of the PU signal can be derived directly from the model assumed for it. As the received power as well 
as the normalized covariance matrix of the PU signal are distinct at each OFDM sub-carrier, PU signal 
detection in this case is executed individually at each sub-carrier, and the results are then combined 
together to form a final decision. In Case B, neither the model of the PU signal nor its distribution is 
known to the detector. The prior knowledge is the frequency band that the PU signal resides. The band 
is assumed to be a continuous segment of sub-carriers. To incorporate the fact that, once a PU signal 
occurs, several sub-carriers in a row are interfered simultaneously, the detector makes its decision by 
jointly considering observations from all possibly interfered sub-carriers. In Case C, no prior knowledge 
of the PU signal is available. Thus, the detection is involved with a search of possibly interfered band. 
The proposed detector is able to detect an abrupt power change when tracing along sub-carriers. 

The organization of this paper is summarized as follows. In Section 11, the signal model of a cognitive 
OFDM system interfered by a PU signal is derived. Three cases concerning the PU signal prior knowledge 
are also described. In Section m, the designs of PU signal detectors are carried out with PU signal prior 
information stated in Section n. Simulation results of the proposed detection algorithms are given in 
Section IV. Finally, we conclude this paper in Section V. 

II. PU Signal Detection in a Cognitive OFDM System 

Consider a wideband cognitive OFDM system with Q sub-carriers. The binary data stream generated 
from the source is encoded and interleaved, and then subdivided into groups of B bits used to generate 
blocks of Q symbols, where each symbol assumes one of L possible values with B = Q log2 L. It is 
assumed that (log2 L)-ary phase shift keying (PSK) modulation is employed. We denote the constellation 
points corresponding to the n-ih block of Q symbols by «S(n) = {S'o(n), Si{n), • • • , S'Q_i(n)}. The n-th 
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TABLE I 

Prior knowledge of the PU signal 





Received Power 


Signal Model 


Gaussian Distribution 


Frequency Band 


Case A 


No 


Yes 


Yes 


Yes 


Case B 


No 


No 


Not necessarily 


Yes 


Case C 


No 


No 


Not necessarily 


No 



OFDM symbol is generated by feeding «S(n) into a Q-point inverse discrete Fourier transform (IDFT) 
and pre-appending the output with cycUc prefix (CP). The resultant signal is up-converted to the carrier 

frequency, and then transmitted over a wireless fading channel. 

At the receiver, after the frequency down conversion and the CP removal, the output signal is passed 
through a Q-point discrete Fourier transform (DFT). In the presence of a PU signal, the DFT output 
corresponding to the n-th OFDM symbol is given by 

Yg{n) = Hg{n)-Sq{n) + Ig{n) + Wg{n), < q < Q - 1, (1) 

where Hq{n) is the frequency response of the channel at sub-carrier q experienced by the n-th OFDM 
symbol, and {Iq{n)} and {Wq(n)} are the contributions resulting from the PU signal and additive white 
Gaussian noise (AWGN), respectively. 

Suppose that a PU signal occupies the frequency band extending from the go-th to the gi-th sub-carriers 
of the OFDM system. If the information of the PU signal frequency band, i.e. qa and qi, is known to 
the detector, the detection algorithm decides whether the signal {Iq{n)} is present in (1) based on the 
observation {Yq{n) : < n < N — 1, qq < q < qi}, where N is the observation length at each sub- 
carrier, and, if any, the prior knowledge of the PU signal. When go and qi are not known, the observation 
{Yq(n)} needs to be extended to all sub-carriers < q < Q — 1. 

TABLE I lists three cases regarding the amount of prior knowledge about the PU signal, including 
the received power, the signal model, probability distribution, and the frequency band it resides. In all 
three cases, the received power of the PU signal is unknown. In Case A, it is assumed the model of 
PU signal is known. Examples that the PU signal characteristic is known to the detector can be found 
in, e.g. [22,26,11]. We assume the sub-carrier indices [go)9i] occupied by the PU signal are known, 
the stochastic process {Iq{n)} observed at each sub-carrier qo < q < qi is Gaussian, and the N x N 
normalized co variance matrices Cg's of the random signal {Iqin)}^^^ at go < g < (?i can be obtained 
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from the PU signal model. The normalization factor to obtain Cg is the PU signal received power at that 
sub-carrier, and Cg has diagonal components equal to one. In Case B, the assumptions of known PU 
signal model and Gaussian distribution are removed. It will be clear this case serves as an intermediate 
stage for developing the detector in Case C, where no prior PU signal knowledge is available. 

III. Design of PU Signal Detector 
A. Case A: Known PU Signal Model, Probability Distribution, and Frequency Band 

Two PU signal models, i.e. a sum of tonal signals and an auto-regressive (AR) stochastic process, 
are used as examples for the detection problem. The PU signal {Iq{n)}^^Q seen at the q-th OFDM 
sub-carrier has a covariance matrix Pi{q)Cq, where Pi{q) is the unknown received power of the PU 
signal at sub-carrier q, and is the normalized covariance matrix with unit diagonal elements. 

1) Tonal PU Signal: Here we model the PU signal as the sum of a number of complex sinusoids. 
Examples include the worldwide interoperability for microwave access (WiMAX) and wireless local area 
network (WLAN) systems, which also employ OFDM technologies. With this model, the received PU 
signal is i{t) = Yli^-oo^l(^ — ITiY, where Tj is the symbol duration, and ii{t) is the signal containing 
the l-th symbol. We have 

K-l 

ii{t) = J2 5?{d,,fe(t)e^(2"^^*+<^)}, 0<t<Ti, (2) 

fc=0 

where K, fi and cf) are the number of complex sinusoids, the carrier frequency, and the random carrier 
phase, respectively, 3fJ{-} denotes the real part, and is the complex baseband signal of the k-Ha 

sinusoid. We have 

di,k{t) = Ci,k ■ Xi^ke'^^'^'/^' , 0<t<Ti, 

where X; ^ is the PSK modulated data of the k-th sinusoid at the l-th symbol, and is the channel 
fading coefficient of the k-th sub-carrier when symbol / of the PU signal is received. We assume, for 
each particular random variables Q^k, k = 0,1, - ■ ■ ,K — 1, are identically distributed. 

Let ri = [Tj/TsJ, where Tg is the symbol duration of the cognitive OFDM, [x\ is the largest integer 
no greater than x, and Pk,q = [(/? — fs + k/Ti)Ts — q]/Q. It is shown in Appendix I-A that, the (n, m)-th 



'Here i{t) is the time-domain PU signal, wliereas {Iq{n)} given in (1) is a frequency-domain signal. 
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element of the normalized covariance matrix Cq{n,m) of {Iq{n)}^^Q is given by 



Cq{n,m) 



E 



j27r(n-m)T,(/i+it/ri) Sm^(7r/3fc,gQ) 



n-m \ k=o - I 1^1 
1 : ' \n-m\<ri-l, 



0, otherwise. 
2) AR PU Signal: We consider the time-domain discrete-time PU signal {ip} at the output of the 
sampler following the frequency down-converter. Suppose that {ip} can be modeled as an r-th order AR 
random process of 

r 

ip = ^ ^ ^j^p—j "I" ^pi (^) 

where {cp} is a white Gaussian random sequence with variance u^, and parameters {<pj}^=i are obtained 
when {ip} has unit power. The contribution of the PU signal at the g-th DFT output for N OFDM 
symbols is Iq = [/^(O), 7^(1), ■ ■ ■ , Iq{N - given by 

I9 = F,i, (5) 
where = diagjf^, fq,-- - , fg} with = [e-i'^^i-^/Q g-^'^'^^VQ . . . g-^'^^^W-i)/^], and i = [io, ii, 

N times 

■■■ )^QAf-i]^- It is readily seen that Iq forms a Gaussian random process as {cp} is modeled to be 
Gaussian. In Appendix I-B, we show how the normalized covariance matrix Cq of Iq given in (5) can 
be computed. 

It is seen that, for the two PU signal models presented above, the normahzed covariance matrix Cq at 
each sub-carrier q G [goj^i] are distinct. Moreover, the PU signal at various sub-carriers have different 
received powers. It will be shown later, cf. (12), these unknown received powers need to be estimated. 
Thus, we conduct the PU signal detection individually at each sub-carrier, and the final decision is made 
by combining the individual decisions at sub-carriers. If a joint detection of all sub-carriers is performed, 
it is required to estimate all unknown PU signal powers jointly. To design a detector based on the 
Neyman-Pearson philosophy, the detection threshold 7^ at sub-carrier q is determined by a given overall 
(i.e., combined from all sub-carriers) false alarm probability PpA = « such that the overall detection 
probabiUty Pq is maximized. Let the decisions made at individual sub-carriers be combined by an OR 
operation, i.e., the detector decides Hi if any of the sub-carriers declares an PU signal is present. For 
both PpA and Pd, we have 

Ps = l- n (l-^s(9)), SG{FA,D}, (6) 

9e[go,gi] 
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where Ps{q) is the detection or false alarm probability at sub-carrier q. Letting -Pfa(q) equal for all q's, 
we obtain the false alarm constraint at each sub-carrier as 

Pfa(9) = 1-(1-«)^/^"", (7) 

where we define the bandwidth of the PU signal as Bpu = qi — q^ + 1. 

Suppose that the detection is performed when the cognitive OFDM system is not transmitting signals. 
The hypothesis testing at the q-ih sub-carrier is 

n = 0,1- ■ ■ ,N — 1, (8) 

Hi: Yq{n)=Iq{n) + Wq{n), 
where {Wq{n)} is complex white Gaussian noise independent of {/^(n)} with distribution CAA(0, cr^I), 
and Ho and Hi represent that the PU signal is off and on, respectively. The PU signal {Iq{n)}^~Q is given 
in (39) and (5), respectively, when it is modeled as a sum of tonal signals and an AR random process. 
Due to the absence of the OFDM signal, {Yq{n)} is a Gaussian random process in either hypothesis. 

Since the detection algorithm proposed for this case is done individually at sub-carriers, we omit q in 
Pi{q) for notational simplicity. The likelihood ratio associated with (8) is 

where = [Yq{0), Yq{l), ■■■ , Yq{N - 1)]^, p{Yq; Pi, Hi) is the probabiUty density function (PDF) of 
Yq under Hi given as 

and p{Yq;Ho) is the PDF under Hq obtained by setting P/ in (10) to zero. Using matrix inversion 
lemma, we have the test statistic InL(Yg) expressed by 

lnL(Y,) = u^PiY\Cq{PiCq + a^I)-'Y, - Indet {PiCq + cr^l) + Incr^. (11) 

It is seen that the unknown P/ in (P/Cg + cj^I)^^ cannot be decoupled from the observation Y^. Thus, 
uniformly most powerful (UMP) test does not exist. Consequently, a generalized likelihood ratio test 
(GLRT) is employed, where P/ in (11) is replaced with its maximum likelihood (ML) estimate P/. 

Let Cq be eigen-decomposed as Cq = VgAgVj, where Vq = [v^^o v^^i • ■ ■ Vq,jv-i] and = 
diag(Ag,o, \i, ■ ■ ■ , \n-i)- Hence, 

det(P7C, + a^,!) = Yl {PjXq^i + a^) and {PjCq + a^^^y' = Vg(PrA, + a^w^y'^l 

i=0 
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The ML estimate of Pj is obtained by substituting the above two relations into p{Yg; Pj,Hi) and finding 
it maximum. Moreover, we should also note that Pj is non-negative. Thus, we have 

Pi = max (o, arg min ^ f ln(PA,,, + a^) + J"""''^'', ^ ^ . (12) 

The general solution of the optimization problem in (12) is unknown, and numerical methods are normally 
required. Even if we can solve (12), the statistic distribution of the detector in (11) is intractable, which 
yields threshold determination of the detector very difficult. On the other hand, under Ho, the random 
variable governing the statistics of Pj is zero half of the time and Gaussian for the other half^. This is in 
contrast to the usual Gaussian asymptotic statistics of an ML estimate. Thus, the asymptotic chi-squared 
distribution of GLRT when N ^ oo does not hold for (11) [28]. 

Due to the difficulties encountered by GLRT stated in the previous paragraph, we resort to a locally 
most powerful (LMP) detector [29,30]. We rewrite (8) as ~ CM{0, PjCg + a^I) with 

no-Pi = versus Hi : Pi > 0. 



The LMP detector, given by 

d\np{Yg;Pi) 



= -a^ tr(C,) + a^^YjC.Y,, (13) 

p,=o 



dPi 

is optimal when Pi is small. Thus, the detector is 

T^(Y,) =YtC,Y,|'7g (14) 

as the remaining part of (13) can be absorbed into the threshold, where the subscript of T^( ) indicates 
it is for Case A. It is seen that LMP has an advantage that no estimate for Pi is needed. Moreover, it 
is almost optimal in the low signal-to-noise ratio (SNR) region for which signal detection is inherently 
a difficult problem. For large departure of Pi from 0, there is no guarantee of LMP's optimality, and 
a GLRT would perform better. However, due to the large SNR, the LMP detector can generally satisfy 
the system requirement with the advantage of lower complexity. An interesting interpretation of LMP 
detectors as covariance sequence correlators can be found in [30, pp. 80]. 



^It is known that, if the probability density function p(x; ff) of the observation x satisfies some "regularity" conditions, then 
the ML estimate of an unknown parameter is unbiased and asymptotically Gaussian (see e.g. [27, Theorem 7.1]). Thus, when 
A'^ is large, the ML estimate in the second argument of max(-, •) in (12), denoted by Pi, is Gaussian. Since Pi — under TCo, 
Pi has zero mean and is larger and smaller than with equal probabilities. It follows that, when N is large. Pi is zero half of 
the time and Gaussian for the other half. 
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Denote by TA(Yy)|w. the shorthand for TA{Yg) under Hi- Let Wg = [Wg{0), Wg{l), ■■■ , Wg{N - 
1)]^. Under Ho, elements of Yg = Wg are independent, and 

N-l 

TA{Yg)U = WtC.W, = ^ A,,|vt^,W,|2 

i=0 

is a weighted sum of independent chi-squared random variables. No general closed form is known for 
its distribution [30, pp. 74-75]. Thus, we look for its asymptotic distribution. We have 

dhip{Yg{ny,Pi) 



dlnp{Yg; Pi) 



E 



P,=o —0 dPi 



, under Ho, (15) 

Pi=0 



dPi 

which by central limit theorem becomes Gaussian. Thus, 

TA{Yg)\no ^Af{a^wN,a^tT{Cl)) , (16) 

where ~" indicates the sense of asymptote, and the formula of 

E{xt Axx+Bx} = tr(AC) tr(BC) + tr(ACBC) (17) 

for X --^ CM{0, C) and Hermitian matrices A and B [31] is employed. Under Hi, due to the PU signal, 
elements of may not be independent. This makes the distribution of T^(Yg)|?^^ difficult to analyze. 
In Appendix 11, the asymptotic distribution of TA{Yg)\-]-i^ is examined using the central limit theorem of 
an m-dependent sequence. 

To determine the threshold 7^ in (14), we represent the cumulative distribution function (CDF) of 
TAiYg)\n, as CDF(x). By (7), jg is given by 

7, = CDF-i((l-a)i/^-), 

with an overall false alarm probabihty a. If TV is large enough such that the asymptotic distribution (16) 
holds, 7g can be further written as 



7, = a^iV + a2^^tr(C2) • Q-\l - (1 - a)i/^™), 

where Q{x) is the Gaussian right-tail probability. On the other hand, if the asymptotic distribution of 
TA{Yq)\T-ig is not valid, histograms of TAiYq)]-]-^^ can be obtained by simulations to get an estimate of 
CDF(x). We can simply produce the histogram with a^r = 1. For any particular a^r, the corresponding 
histogram can be easily mapped from that of a^r = 1. 
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PU signal detection ^ 


1 ^ 


i i 




Channel estimation 






1 




Payload 











r = 



Fig. 1. The signalling of the cognitive OFDM system for the detection algorithm developed in Case B. 

B. Case B: Known PU Signal Frequency Band 

In this case, we employ the fact that the PU signal, if present, appears simultaneously at the sub- 
carriers from qq to qi. The algorithm developed for this case works for PU signal with the bandwidth 
^pu = (Zi — ^Zo + 1 > 2. When BpTj = 1, the PU signal can be detected by first estimating its received 
power and then employing an energy detector [4]-[10]. 

Unlike in Case A that the detection is performed when the cognitive OFDM is not transmitting signals, 
the detection method presented for Case B can work when the CU signal is present. The OFDM system 
signalling for the detection algorithm in Case B is illustrated in Fig. 1. Before initiating (t = 0), the 
system performs a PU signal test at the suspect sub-carriers. If the PU signal is present, the cognitive 
system does not send any signal over these sub-carriers. On the contrary, if the PU signal is absent, 
channel estimation of the cognitive system is carried out, and the payloads are then transmitted over 
them. If necessary, during payload transmission, PU signal testing may be executed periodically to ensure 
a quick response to the appearance of the primary network. As shown in Fig. 1, PU signal detection, 
channel estimation and payload transmission are repeated over and over again (adding PU signal test 
during payload transmission, if necessary) until the presence of PU signal is detected. Once PU signal is 
detected either at the system initialization or in the middle of a normal operation, the OFDM system stops 
transmitting signals over sub-carriers q € [qo,qi]- That is, channel estimation and payload transmission 
are suspended, while the PU signal detection is still performed periodically. 

The following situation may arise when PU signal monitoring is done concurrently with the transmission 
of the cognitive system. When a missed detection occurs, the channel estimation will be executed under 
the presence of the PU signal, resulting in a poor estimation accuracy. These inaccurate channel estimates 
are subsequently adopted in the next PU signal detection (see (19)-(21) below), which is expected to 
deteriorate the detection performance. To escape from the vicious cycle, it is required that PU signal 
detection is carried out in the absence of the CU signal. This can be done by enforcing a "silent period" 



October 6, 2008 



DRAFT 



11 

every a given time interval, during which the CU should suspend the transmission. Although the insertion 
of silent periods enables the CU to escape from the vicious cycle, a side effect occurs that the efficiency 
of CU is reduced. A discussion of how often and for how long the cognitive system should remain silent 
in a given transmission interval is an important research topic, e.g. [32]-[35]. A trade-off should be made 
between two opposing issues of efficiency and integrity of sensing results to produce a desirable balance 
defined by a suitable objective function. For example, in [35], an MAC-layer sensing period adaptation 
algorithm is designed to maximize the discovery of spectrum opportunities. 

As described in the previous paragraph, PU signal detection may be executed when the cognitive 
OFDM is either on or off. The received signal at the g-th sub-carrier Yq{n) is given by (1) with Hq{n) 
and Sq{n) set to zero when the OFDM system is not transmitting. For < q < qi, build an 
observation Z{q) from {Xq{''^)}n=Q ^^^^ that the PU signal detection is based on the observation along 
the frequency domain, i.e. Z = [Z{qQ), Z{qQ -|- 1), • • • , Z{qi)]'^. We choose 

Z{q) = ^YlYq, qo<q<qi, (18) 

because Z{q) is the periodogram of the received signal at the g-th sub-carrier averaged over N OFDM 
symbols. It is well known that the periodogram is an estimate of the true spectrum of a signal. Another 
interpretation of (18) is that the normaUzed co variance matrix Cq in the LMP detector of (14) is replaced 
with the identity matrix due to the unavailabiUty of it. That is, elements of are regarded as uncorrelated. 
Expanding (18), we obtain 
1 

n=0 

2R{Hqin)Sq{n)rq{n)} + 2U{Hq{n)Sqin)W;{n)} + 2»{/,(n)W^;(n)}) . (19) 

Depending on whether the OFDM system is transmitting or not, Hq(n) is either known from chaimel 
estimation or equal to zero. We define 

1 

n=0 

where Hq{n) and a^^ are the estimates of Hq{n) and a^^, respectively. A new observation is built as 

Z{q) = Z{q) - m{q), qo < q < qi, (21) 

which corresponds to subtracting the first and third terms inside the brackets of (19). 

Let Z = [Z{qo), Z{qo -|- 1), • • • , Z{qi)]'^ be the observation of the spectrum sensing problem. It is seen 
that each component of Z is involved with a number of contributions from the PU signal, CU signal. 
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AWGN, estimation errors of Hg{n) and cr^, and so on. Thus, it is hard to make a precise statistical 
description for Z, leading to the difficulty in formulating the corresponding hypothesis test. To alleviate 
the problem, here we adopt the model investigated in [36], where there exists an interfering signal lying 
in an arbitrary unknown subspace of the observation space. Specifically, when the PU signal is present, 
the observation vector Z is formulated as 



I + R + N, 



(22) 
(23) 



w-i 



Af-l 



N-l 



where 

I = ^ E Yl l^.o+i(n)P, ■■■,Y1 (24) 

Ln=0 n=0 n=0 

is the contribution of the PU signal, R is the component due to an unknown interference, and N is the 
white noise. The PU signal I resides in an (r + 1) -dimensional subspace spanned by the columns of the 
known Bpu x (r + 1) matrix H, given as 

/ ho{0) hi{0) ■■■ hr{0) \ 

hoil) /ll(l) ••• hr{l) 



H= [ho,hi, 



(25) 



\ho{qi-qo) hi{qi - qo) ■■■ hr{qi-qo) J 
and has an unknown gain vector //. That is, the powers of the PU signal across sub-carriers q G [qo, qi\ 
is modeled as a linear combination of vectors {hj:0<i<r}. The white noise N is a Bpu -dimensional 
Gaussian random vector modeled as tr^N, where the scalar cr^ is unknown, and the covariance matrix 
of N is the identity matrix. Finally, the vector R = accounts for the effects that are ignored by I 
and N in (22); both the matrix U, whose columns constitute the subspace of R, and the gain vector ^ 
are unknown. 

It is argued in [36] that we require to robustly choose the unknown matrix U in the formulation of a 
hypothesis test such that an adequate level of protection to false alarm as well as a sufficient detection 
sensitivity to the signal are both maintained. Let G be a Spy x (-Bpu — r — 1) matrix whose columns span 
the orthogonal complement of the space generated by H, i.e. G = H-*-. A miiumax-based reasoning 
in [36, Appendix A] leads to explicit and different choices for the unknown subspace U in the two 
hypotheses. That is, U = G for TYq, and U is the zero matrix for Tli. Thus, the hypothesis test is 

Hq: Z = GV' + a2N, 
Hi- Z = H/i + (jf N, H/i )^ 0, 



(26) 
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where H/x )^ means that all elements in H/i are non-negative, and /x, ctq and a\ are all unknown. 
To perform PU signal detection, GLRT of 

p(Z;^,c72,7^o) 

is employed, where /t, ^, and cr? are ML estimates of /i, ^ and cr?, respectively. Although, when T^i is 
true, /i and erf are both parameterized by the PU signal; however, joint estimate of these two unknowns 
results in a complex detector structure. Thus, in spite of their dependence, /i and a\ are estimated 
separately. The specification of the PU signal subspace H is only an approximation, and the performance 
of the detector depends on whether the linear subspace spanned by {hj} gives good description of the 
signal class. In the case where the correlation matrix of the PU signal is known, H could be selected 
as orthogonal eigenvectors of the correlation matrix, and the subspace dimension r + 1 can be chosen 
based on some information measures [37], e.g. the Akaike information criterion (AIC) [38] and the 
minimum description length (MDL) [39]. When the correlation matrix of the PU signal is unknown, 
the above information measures and their variations, e.g. [40], can assess the discrepancy between the 
true and approximating models, which serve as useful tools in solving the model selection problem. We 
suppose that the set of vectors {hj} used to model the PU signal chaimel selectivity is suitably chosen. 
Consequently, H/i :>= holds most of the time, and the one-sided test of H/i in (26) does not bring much 
trouble. 

It is shown in [36] that the hkehhood ratio in (27) for the robust hypothesis test leads to the matched 
subspace filter, given by 

T (7\ = ^PU-r-1 Z^H(H^H)-iH^Z 

' r + l Z^(I-H(HrH)-iHr)Z' ^ ^ 

It is known that, for the test in (26), when G is set as the zero matrix, GLRT yields the matched 

subspace filter [41]. Equation (28) demonstrates, even in the presence of unknown G^*, the matched 

subspace detector is optimal, meaning that it is robust to the interference whose subspace is unknown. 

Under Ho, the detector is distributed as Tb(Z) ~ Fj.+i^Bn]-r-i, where Fa^^ denotes an F distribution 

with a numerator degrees of freedom and b denominator degrees of freedom. Given threshold 7, the false 

alarm probability is given by 

with QFa,Ax) the right-tail probabiUty of F^ ;, evaluated at x. If {hj} is able to model the PU signal I 
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well, we have Tb{Z) ~ F^_^_^^^^^_^_^{X) under 7ii, where 

, 91 N-1 

^ = -T.T.\w\'^ (29) 



<7T 

q=9o n=0 



and denotes a noncentral F distribution with a numerator degrees of freedom, b denominator 

degrees of freedom and non-centraUty parameter A. Thus, the detection probabiUty with threshold 7 is 



where Qp' ^{x){x) is the right-tail probability of F^j^{X) evaluated at x. 

C. Case C: No Prior Knowledge of PU Signal 

In this case, the information of possibly interfered frequency band is unknown. Consequently, the 
detection algorithm should be involved with a search of the interfered band. Intuitively, given the 
observation Zo:Q-i = [Z{0), Z{1), ■ ■ ■ ,Z{Q — l)]'^, this search is based on the powers at all sub- 
carriers, and, if the cognitive OFDM system is transmitting signals, a sub-carrier with larger frequency 
response magnitude \Hq{n)\ tends to be judged as the PU signal is present. To avoid this problem, the 
search of interfered band is executed when the cognitive system is not transmitting signals. 

The hypothesis testing associated with Case C is a detection of abrupt changes [42], given as 

Ho: Z{q)^Uo{q), qe[0,Q-l], 

Uo{q), q G [0, go - 1] U [^i + hQ- 1], (30) 

Ui{q), q G [qo,qi], 

where {Uo{q)}q are white Gaussian with variance Uq, and {C^i(g)}g=go independent Gaussian with 
the mean vector modeled by {hj} and variance af. All of qo, qi, (Tq, af and the weighting factors of 
{hi} are unknown. 
We can obtain that the GLRT corresponding to (30) is 

max — — (31) 



Hi: Z{q) 



«o,ai (0-2 ) (Q-ai +ao-l)/2 (o-2 ) (a, -ao+l)/2 



where 



2 

and afi^^ 



{Q-ai+ao- 1) (Zo.^^_^Zo:ao-l + '^a. + l-.Q-l^ai + l-.Q-l), 
(01 -aQ + 1)-^ZT^ (I - Hag;aj(HT^ Hag:aJ-lHT^ )Za 



denote estimate of ctq under T-Lq, estimate of ctq under Wi, and estimate of a\ under H\, respectively, 
and H-ao-.ai is given in (25) with go and q\ replaced by ao and a\, respectively. Defining /(ao, ai) as the 
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target to be maximized in (31), we consider the false alarm probability for the detector with threshold 
7, i.e., 

PpA = Prob <^ max /(ao, ai) > j;Ho} , 

= l-Prob{/(ao,ai) <7,V[ao,ai] C [0, Q - 1]; T^o} • 

Since the random variables governing /(ao,ai) for different choices of ao and ai are not necessarily 
independent, the determination of PpA and hence the detector threshold for a given PpA = a becomes 
intractable. 

To conquer this problem, the PU signal detection is decomposed into two steps. In the first step, we 
search for PU signal's frequency band by Zo:Q-i to get estimates % and qi. In the second step, we 
assume go and qi obtained in the first step are correct, and we can consequently perform PU signal 
detection in the same way as that proposed for Case B, where PU signal frequency band is known. In 
specific, the first step solves the optimization problem of (31). As the numerator is not a function of ao 
and ai, the optimization is equivalent to miiumizing the denominator, i.e. 

{qo,q\) = arg min (ag, )W-«^+«o-i)/2(a2 (32) 

(Oo,Ol) ' ' 

However, this problem is complex because there are about possible trials for combinations of ao 
and ai. 

To reduce the computational load, we can simpUfy the optimization in (32) to one that minimizes the 
least square (LS) error between the observations Z{qys and the estimated PU signal power, i.e., 

(go,9i) = arg min ^ Z{qf + ^ I Z(g) - ^ /xj/tj(g - ao) j . (33) 

Equation (33) is interpreted as follows. Analogously to the discussion in Case B, the vector Z^^ig^ = 
[^(go), -2^(90 + 1)) ■ ■ ■ 5 contains a signal lying in the (r + l)-dimensional vector space spaimed 

by columns of Hq;,:^^ and having an unknown gain vector /i. Consider the second term of the target 
function at the right-hand-side of (33). Supposing the PU signal resides at sub-carriers q G [ao,ai], we 
find the LS estimate of the gain vector /t = [/io,/ii, • • • , Ar]^, i-e. A = ^ao-.aj^a^-ay^^a^-.a^'^a^.o.^, 
and compute the LS error between T^a^-a^ and Ha^^a^/t. On the other hand, in the first term of the target 
function, only the sum of is taken into account since sub-carriers g G [0, Q — 1] \ [oq, ai] are free 

from the PU signal. 

To compare the ML estimates of go and q\ and the suboptimal ones, the target function of (33) is 
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written as the vector form 

to facilitate examining its relation to the target function of (32), where the third term of (34) is the LS 
error yielded by the second term in the target function of (33). 

The solution of (33) can be found by the technique of dynamic programming (DP) [43,44] as follows. 
Define ^ 
5o{a,b):= Z{qf and 5i(a,6):= J] [ Z(g) - ^ - a) ) , 

where {fii}l=Q is the LS estimate that minimizes 6i{a,b). Let 

e{l)= min (5o(0, ao - 1) + (5i(ao, 0, ^ < ^ < <5 - 1, (35) 

0<ao<l—r 

where the constraint Qq < I — r guarantees existence of /xj's in 5i{ao,l). The optimization in (33) is 
equivalent to 

min 5o(0, oo - 1) +Si{ao,ai) + do{ai + 1,(5-1) 

(aojOi) 



= min< ( min5o(0, ao - 1) + Si{ao,ai) ] + do{ai + l,Q -1] 
= min e(ai) + So{ai + 1,Q - 1). (36) 

r<ai<Q—l 

Thus, qi can be found by searching for ai that minimizes (36), and go is equal to the value of the 
argument ao in (35) that minimizes e{qi). Note that, the computation of 6i{ao,l)'s in solving (35) can 
be done recursively by sequential LS formulas [27, pp. 242-251], and the DP can reduce the complexity 
of search from the order of to the order of Q. 

Suppose that DP yields correct values of qo and qi. We then employ the detector proposed for Case 
B to decide whether a PU signal is present in the estimated frequency band. The performance analysis 
of the detector is executed as follows. The false alarm probabiUty Pfa of the detector in Case C is 

^ Prob {go = ao, qi = ai; Wo} • Prob {TB{Zaa.,a,) > lao,a,;'Ho} , (37) 

[ao,ai]c[0,Q-l] 

where the first probabiUty is the one that DP yields the result of (go, gi) = (ao, ai) under Ho, and Tb{-) 
is given in (28) with Bpn and H there replaced by ai — ao + 1 and Hao-.ai, respectively. We caimot 
determine the first probability in (37). However, for each DP searching result (go, gi) = (ao, ai), we can 
choose a threshold 7ao,ai for TB(Za„:aJ such that the second probability in (37) is a constant. In this 
case, the false alarm probability Pfa becomes 

Pfa = Prob {rB(Z ) > 7ao,ai;'Ho} ■ 
October 6, 2008 DRAFT 



17 



Given a constraint of PpA = we choose 

7oo,oi = [ao, oi] C [0, Q - 1]. 

Detection occurs when, under Tii, DP returns a correct result and the decision statistic is greater than 
the threshold. Let p denote the probability that the frequency band search returns a correct result. When 
the PU signal is well modeled by {hj}, the detection probability Pq is equal to 

p-Prob{Ts(Z,„.,J > Q-pl^^^^.,„Ja);m] 

where A is given in (29). 

IV. Simulation Results 

Throughout the simulations, the tonal model presented in Paragraph III-A. 1 is adopted for the PU 
signal. The parameters of the PU and cognitive OFDM systems are Tg = 312.5 ns, Ti = 26.6 ;us, 
fs = 3.1 GHz, fi = 3.36 GHz, and Q = 128. The number of complex sinusoids K contained in the PU 
signal is adjusted according to the PU signal bandwidth Bpu. 

In Fig. 2, the receiver operating characteristic (ROC) of several detectors are shown to illustrate the 
performance of the detector proposed in Case A. The probabilities of miss and false alarm are shown in 
the vertical and horizontal axes, respectively, where the detection is performed at a single sub-carrier. The 
overall detection performance considering all sub-carriers can be obtained by (6). The observation length 
N is set to 80 OFDM symbols. During the detection, the cognitive system is not transmitting signals. 
The curves in the figure are divided into two groups for the power ratio of the PU signal and AWGN 
as and —2 dB. Within each group, there are four curves. The three solid ones from top to bottom are 
simulation result of energy detector (test statistic YqYq), simulation result of LMP detector, and analytical 
result of estimator-correlator^, respectively; the dashed line is the analytical result of LMP yielded by 
Gaussian approximation, i.e. r^(Yq)|'^(, and r^(Yg)|'^^ are distributed as (16) and (48), respectively. 
Consistently with our intuition, the estimator-correlator has the best performance due to its full knowledge 
of the observation's PDF, and the energy detector is the worst since the correlation in the PU signal is 
not exploited. In getting the distribution of (16), the central limit theorem for independent and identically 



^The estimator-correlator is derived from the likelihood ratio test with known received power of the PU signal. The estimator- 
correlator and its performance can be found in [29, pp.142]. 
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10 




P , false alarm probability 



Fig. 2. ROC curves for energy detector, LMP detector and estimator-correlator performed at a single sub-carrier; Q — 128, 
A'^ — 80, quiet cognitive system. 

distributed random variables is adopted, while a central limit theorem for m-dependent random variables 
is used to arrive at the distribution in (48). The significant discrepancy in the simulated and analytical 
results of LMP detector demonstrates the assumptions in achieving the Gaussian approximations, in 
particular for (48), are not valid under the simulation environments. 

In Figs. 3 and 4, the simulated ROC curves of the detector proposed in Case B, i.e. (28), are plotted for 
the environments that the PU signal experiences a channel of eight multipaths with uniform power delay 
profile and IEEE 802.15.3a ultrawide band (UWB) CM3 model, respectively. In either case, a number 
of channel realizations are run to obtain an averaged performance. The function /ij(n) in (25) is set as 
71*. That is, an r-th order polynomial is used to model PU signal powers across the sub-carriers. The 
observation length is 70 OFDM symbols. During detection, the cognitive system is quiet. The average 
power ratio of the received PU signal and AWGN over the affected sub-carriers is controlled to be dB. 
The bandwidth of the PU signal Bpu and the order r of the polynomial used to model the PU signal 
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10 10 10 10 



Pp^, false alarm probability 

Fig. 3. ROC curves for detector proposed in Case B when the PU signal experiences a multipath channel with path number 
equal to eight and uniform power delay profile; Q = 128, = 70, quiet cognitive system, average PU signal to AWGN power 
ratio dB. 

powers are indicated on each curve. It is shown that, given the same values of Bpu and r, the detector 
performance shown in Fig. 3 is better than that in Fig. 4, indicating a severe frequency-selective channel 
of the PU signal deteriorates the performance. This is because the polynomial fails to model the PU 
signal powers in a hostile channel. Another observation is that, under the same channel type, detection 
performance improves when i?pu increases; whereas increasing the polynomial order r is not necessarily 
helpful for detection. This can be explained as follows. The numerator and denominator of Tb(Z) can be 
seen as estimates of the PU signal power and noise power, respectively [41]. The dimensions of the signal 
and noise subspaces are r + 1 and Spu — r — 1, respectively. If r is too large, the signal is overestimated, 
which increases the false alarm probability; on the other hand, if r is too small, the probability of miss is 
increased. As mentioned in Section III-B, information measures can adopted to choose a suitable value 
of r. 
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Fig. 4. ROC curves for the detector proposed in Case B wlien thie PU signal experiences an IEEE 802.15.3a UWB CM3 
channel; Q = 128, N = 70, quiet cognitive system, average PU signal to AWGN power ratio dB. 

In Fig. 5, the performance of the detector proposed in Case B is demonstrated when the cognitive 
system is transmitting signals. An r-th order polynomial is used to model the PU signal power. The 
power ratio of the CU signal to AWGN is set to 8 dB at every sub-carrier. The observation length is 70 
OFDM symbols. In Fig. 5(a), the average PU signal to noise power ratio versus is plotted when the 
PU signal experiences a multipath channel with path number equal to eight and uniform power delay 
profile. The bandwidth i3pu, P^a, and polynomial order r are indicated on each curve. Fig. 5(b) shows the 
same information as Fig. 5(a) with the channel experienced by the PU signal being an IEEE 802.15.3a 
UWB CM3 channel. 

In Fig. 6, the results of interfered band search based on (33) of Case C are shown. We try different 
Spu and observation length N with a constant PU signal power over the sub-carriers, i.e. a frequency-flat 
fading channel. The dimension of signal subspace is set to 1, and ho is an all-ones vector. A search is 
regarded to be a hit if |go — go| < 1 and jgi — gi| < 1. It is shown that, with an observation length 
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SNR (dB) 



(b) 

Fig. 5. Detection probability versus the power ratio of the PU signal to AWGN for the detector proposed in Case B, when 
the power ratio of the cognitive signal to AWGN is 8 dB, and the channel experienced by the PU signal is (a) a multipath 
channel with path number equal to eight and uniform power delay profile, (b) an IEEE 802.15.3a UWB CM3 channel; Q = 128, 
N = 70. 
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Fig. 6. Hit percentage of the interfered band searching versus the power ratio of PU signal to AWGN for different values of 
Bpu; constant PU signal power at all sub-carriers, Q — 128, quiet cognitive system. 



N = 70, the hit percentage approaches 100% at about SNR 1 dB. It is also seen that the hit percentage 
is irrelevant to the bandwidth Bp\]. 

However, when the PU signal experiences severe frequency-selective fading channel, where there is 
a notch within the interfered band, our simulation results indicate that the performance of band search 
is not as good as that presented in Fig. 6. Examples of channels that result in erroneous estimates are 
shown in Fig. 7, where the responses (sub-carriers 30 to 39) at which the PU signal resides are plotted. 
For each sub-figure, the horizontal and vertical axes are sub-carrier index and squared magnitude of 
channel frequency response, respectively. In Figs. 7(a) and 7(b), there are spectrum notches within the 
band; while, in Fig. 7(c), all the sub-carriers suffer from deep fades. In the cases of Figs. 7(a) and 7(b), 
the estimate [qo,qi] is a subset of the true band [gcgi]; in the case of Fig. 7(c), [go,gi] is not even 
a subset of [gcQi]. Such problem arising in a severe frequency-selective channel can be conquered by 
cooperation among CUs, as channels from the PU signal source to CUs of non-proximity can be regarded 
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Fig. 7. Examples of frequency-selective fading channels that result in erroneous estimates of interfered band. The vertical axis 
is the squared magnitude of the channel frequency response. 



as independent. 

V. Conclusion 

In this paper, the problem of PU signal detection in an OFDM based cognitive system is addressed. We 
categorize the amount of PU signal's prior knowledge into three cases. In each case, a Neyman-Pearson 
detector that exploits the prior information is proposed. 

In Case A, it is assumed that the PU signal model is known, and the PDFs of the received signal 
under both hypotheses are completely available except for the received power. The frequency band that 
the PU signal resides is known as well. Due to the difficulties of finding the received signal power 
estimate and the detector threshold, the use of a GLRT detector is not suggested. In stead, an LMP 
detector is employed due to the advantages of being optimal for weak PU signal and no need to get the 
unknown power estimate. Since the covariance matrix and the received signal power are distinct at every 
sub-carrier, the detection is performed individually at sub-carriers, and the final decision is formed by 
an OR operation of the results of the individual detections. The simulation result of an LMP detector 
is compared with an energy detector (information of PU signal covariance matrix not exploited) and 
estimator-correlator (perfect knowledge of PU signal received power). The performance of LMP detector 
is between the other two. 

In Case B, we assume the only PU signal prior information is its residing sub-carriers. The detector 
proposed for this case exploits the fact that, once the PU signal appears, it interferers a consecutive 
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segment of sub-carriers simultaneously. Measurements obtained at all these sub-carriers are taken as the 
observation. The hypothesis test is an unknown subspace signal detection in an unknown interference and 
a Gaussian noise with unknown variance. The result of [36] is employed to robustly choose the subspace 
where the unknown interference locates, and it is shown a matched subspace detector is the GLRT of 
the hypothesis testing. 

In Case C, no prior knowledge about the PU signal is available. The detection is involved with a 
search of interfered sub-carriers. It is shown that the GLRT has high complexity and is difficult to find 
the output distribution, leading to an undetermined threshold. Thus, the detection is divided into two 
steps. The first step searches for the interfered band using an ML criterion, and the second step decides 
whether a PU signal is present in the estimated interfered band. It is seen that the second step corresponds 
to the problem considered in Case B. The first step is further simplified by employing an LS criterion 
instead of ML, which enables the use of the DP technique to solve the optimization problem. Simulation 
results show the search of interfered band has a very high accuracy if the channel experienced by the 
PU signal is frequency-flat faded. However, if the channel is severely frequency-selective, the estimation 
accuracy degrades. It is beHeved that cooperation among CUs is helpful in conquering this problem. 

Appendix I 

Normalized Covariance Matrices of PU Signal Models 
A. Tonal PU Signal 

Let the cognitive OFDM system have symbol duration Tg, carrier frequency fg, and its zero-th symbol 
start at t = to- We suppose that, at the time the n-th OFDM symbol is received, i.e. nTs + to<t< 
{n + l)Ts + tQ, it is within the span of PU signal's ^-th symbol.'* When detecting the n-th OFDM symbol, 
at the down-converter output of the receiver, the contribution resulting from the PU signal is given by 

K-l 

k=0 

It is then sampled every := Tg/Q seconds, resulting in Q samples of 

ip = E 0,.^/,fee^''-(^W^O^^^e^-^(-'''^), 0<p<Q-l, (38) 

k=0 



Although it is possible that an OIDM symbol crosses the boundary of two PU signal symbols, as T; is in general much 
larger than Tg, the contribution resulting from this situation is negligible. For instance, if Mobile WiMAX and muhiband (MB)- 
OFDM-UWB are the sources of the PU signal and the cognitive OFDM system, respectively, the value of Ti/Ts is as large as 
329. 
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where A/ = fi — fs and 0{n, I, k) = 27r[(nTs + to)(/i + k/Ti) — IfiTi] + (p. The discrete-time signal 
{ip} is passed through a Q-point DFT, which gives 

Q-l 

p=0 

= E 0,,X,,e^-^(-'''^)e^-/^-^^(Q-^) , 0<g<Q-l, (39) 
fc^ sin(7r/3fe,5) 

where ^k,g = (A/ + A;/T,)rd - q/Q. It is seen from (39) that, due to central limit theorem, {Iq{n)}^_Q 
can be approximated by a Gaussian random sequence when the number of complex sinusoids K is large 
enough. 

Let us denote by event A that the n- and m-th OFDM symbols both fall within the span of the Z-th 
PU signal symbol. Clearly, Iq{n) and Iq{m) are zero-mean random variables. Conditioned on event A, 
the correlation of Iq{n) and Iq{m) is 

E{Iq{n)i;{m)\A} = E{|C,,P}e^-^-(--)^-/' ,.2.(n-^).T./r. sin^(^/?MQ) 

where the expectation at the left-hand-side is with respect to symbols X; ^ of the PU signal and the 

fading coefficients Ci,fe- On the contrary, if OFDM symbols n and m fall within the spans of distinct 
symbols of the PU signal, we have 

E{Iq{n)I*g{m)\A} = 0, 

where A is the complement of A. 

Let r) := [Ti/T^J with [xj the largest integer no greater than x. The probability Pr{^} of event A is 
roughly equal to 



Pr{^} = < 



_ n — m I I 

1 , In — m\ < rj — 1, 



0, otherwise, 

where we omit the case the OFDM symbol(s) fall at the border of two PU signal symbols. Thus, the 
(n, m)-th element of the covariance matrix of {Iq{n)}^^Q is given by 

Pr{^}.E{I,(n)/*(m)|^}, 

which leads to the normalized covariance matrix Cq in (3) assuming, for each particular I, 0,it's are 
identically distributed. 
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B. AR PU Signal 

The method to obtain the normalized covariance matrix for AR PU signal basically follows the hne 
presented in [45, pp. 414]. Since Iq has zero mean, the covariance matrix of Iq is given by 

E{I,lt} = FgRiFt, (40) 

where ^ denotes Hermitian transpose, and Ri is the correlation matrix of the time-domain PU signal 
signal i. From (4), it is readily seen that 

/ 1 \ f io \ f io 



(f>r 



<t>r 



\ 



h-1 



(41) 



V 0r • • • ^1 1 / V iQN-1 / \ eQN-1 / 

which has a notational form of 

io:r-l 1 _ ( io:r-l 

ir:QN-l J \ ^r-.QN-l 

where notation definitions can be easily mapped between (41) and (42). Right-multiplying both sides of 
(42) by their Hermitian transposes and then taking expectations, we obtain 



(42) 



\ i/^I 

where Rio = E{io:r— iIq ,.-1} is unknown. Since A is non-singular, we have 



(43) 



Rr 



1 







I A. 

u-n 



(44) 



Rfi = At 

Partition A into the following four blocks: 

A = 1 , (45) 

y A21 A22 

where the dimensions of A21 and A22 are (QN — r) x r and {QN — r) x {QN — r), respectively. 
Plugging (45) into (44), we can write 

1 _ / ^iL. + ^"'aLA21 .^-^At^ A2 



R." 



^22 



V 2AI2A21 



(46) 
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It is known that the inverse of a nonsingular ToepUtz matrix is persymmetric (i.e., symmetric about 
the northeast-southwest diagonal) [46]. Since Ri is ToepUtz, Rio^_i can be determined by comparing 
z/~^A2iA2i and z/~^A22A22 in the northwest and southeast blocks of the block matrix in (46). Then, 
the right-hand-side of (40) can be written as 

F,A-M ^'"^ ^ ](At)-iFt. (47) 
\ j 

Note that since the parameters are obtained when the PU signal {ip] has unit power, the matrix 

given in (47) is the normalized co variance matrix of Iq. 

Appendix II 

Asymptotic Distribution of LMP Detector Output under Hi 

To consider the asymptotic distribution of the LMP detector output under Hi, we employ a central 
hmit theorem for an m-dependent sequence. The following definition and theorem are helpful. 

Definition 1: [47, pp. 69] A sequence of random variables, Xi, X2, ■ ■ ■ , is said to be m-dependent if 
for every integer, s > 1, the sets of random variables {Xi, • • • ,Xs} and {Xjn+s+i^ Xm+s+2^ ■ • ■ } are 
independent. ■ 

Theorem 1: [47, pp. 70] Let Xi,X2, - ■ ■ , be a stationary m-dependent sequence with finite variance 
and let 5„ = EILi Then [Sn ~ E(5„)]/ Vvar(5„) converges in distribution to J\f{0, 1). ■ 

In the following, we will show that the LMP detector output YqCqYq under Hi is a sum of random 
variables in an m-dependent sequence. We have 

N-l N-l N-l N-l 

YlCqYq =J2Y1 Yq{zrCq{i,j)Yq{j) = ^ Yqii)* ^ Cq{z, j)Yq{j) , 

where Cq{i,j) is the (i, j)-th entry of Cq. Define a random variable Xi as 

N-l min(t+r;— l,Af— 1) 

X, = Yq{i)* Yl Cq{i,j)Yq{j) = Yq{i)* ^ Cq{i, j)Yq{j) , 

i=0 j=max(0,j— »7+l) 

where tj satisfies that Cq{i,j) « whenever \i — j\ > rj. In the case of a tonal PU signal model, we have 
rj = [Ti/Ts\. It is clear to see {X^} is an m-dependent sequence, and YqCqYq is a sum of m-dependent 
random variables. By Theorem I, if {Xi} is stationary and each X^ has finite variance, then the LMP 
detector output under Hi, i.e. r^(Yg)|7^j, is asymptotically Gaussian. At this moment, we have 

TA{Yq)\H, ^M{tr{Cq{PqCq + CT^I)}, tr{(C,(P,C, + CT^I))'}), (48) 

where the statistical property of Gaussian in (17) is used. 
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